Dephasing via stochastic absorption: A case study in Aharonov-Bohm oscill ations 



O 

o 

(N 

U 

Oh- 
< 



Colin Benjaminf] and A. M. Jayannavar[j] 

Institute of Physics, Sachivalaya Marg, Bhubaneswar 751 005, Orissa, India 
(Dated: February 1, 2008) 

The Aharonov-Bohm ring has been the mainstay of mesoscopic physics research since its inception. 
In this paper we have dwelt on the problem of dephasing of AB oscillations using a phenomenologi- 
cal model based on stochastic absorption. To calculate the conductance in the presence of inelastic 
scattering we have used the method due to Brouwer and Beenakker. We have shown that conduc- 
tance is symmetric under flux reversal and visibility of AB oscillations decay to zero as a function 
of the incoherence parameter thus signalling dephasing in the system. Some comments are made 
on the relative merits of stochastic absorption with respect to optical potential model, which have 
been used to mimic dephasing. 
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Dephasing is denned as the process by which auantum 
mechanical interference is gradually destroyecoB. De- 
phasing of electronic phase coherence may be caused by 
it's interaction with other quasi particles in the system. 
The notion of intrinsic decoherence and dephasing of a 
particle interacting with it's environment is being investi- 
gated intensively. For, a comparison between theory and 
experiment, it is necessary to know how dephasing affects 
the various quantum phenomena related to transport. In 
our present work we will be interested in modeling de- 
phasing phenomenologically. The Aharonov-Bohm (AB) 
interferometer is one of the best examples for analyzing 
how quantum interference effects are affected by dephas- 
ing. In this interferometer the phase of the electrons pass- 
ing through the arms is modulated by the magnetic flux. 
It is not unlike the Young's double slit experiment apart 
from the presence of the magnetic flux. In the Young's 
double slit interferometer as we know the intensity is 
given by I = |#| 2 = |^i \ 2 + \ip 2 \ 2 +2Re(ipl'if; 2 e l ' i3 ) , the part 
2Re(ipiip 2 e ) represents the interference term. Here 
and ip2 are the complex wave amplitudes across the upper 
and lower arms of the interferometer (this discussion is 
restricted to only two partial amplitudes for simplicity) 
and 4> is phase difference between the two wave ampli- 
tudes. If there is no phase relationship between the waves 
then the intensity will be I = l^il 2 + IV^I 2 - Dephasing 
thus suppresses the interference terms. In other words 
dephasing, leads to the vanishing of the off-diagonal ele- 
ments of the density matrix. 

Aharonov-Bohm effect in mesoscopic rings manifests 
itself as a periodic behavior of conductance as a function 
of magnetic flux piercing the loopS. In an AB interfer- 
ometer the electrons are treated quantum mechanically 
along the ring and they retain their phase memory across 
the entire sample and this gives rise to an operative def- 
inition of the phase coherence length or incoherence 
length. On increasing the temperature becomes less 
than the size of the ring and the phase coherence in the 
system gradually disappears and as a consequence the 
AB effect (the visibility of the AB oscillations) vanishes. 
Typically scales with temperature T in a power law 
form, i.e., L$ — T~ p (p lies in the range 1 to 2). 



There are different ways to model dephasing in such 
mesoscopic systems. A controlled way to introduce de- 
phasing in the AB interferometer is to attach a voltage 
probetl to the sample as in inset of Fig. 1 (Buttiker's 
model). This voltage probe breaks the phase coherence 
by removing electrons from the phase coherent channel 
in the system and subsequently re-injecting them with- 
out any phase relationship. Another way to introduce 
dephasing, is to add a spatially uniform imaginary ppr 
tential (— iVi) (or optical potential) to the Hamiltonianta, 
which removes particles from the phase coherent motion 
in the system. In this case absorbed particles are iden- 
tified as a spectral weight lost in the inelastic channels 
and they are re-injected back into the sample. In the 
method due to Zohta and EzawaQ the total transmission 
is defined after re-injection as the sum of two contribu- 
tions one due to the coherent part and the other due to 
the incoherent part, i.e., T tot — T C( 
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herent part is calculated as Ti nco h = Ti+t 
and Ti are the probabilities for right and left transmis- 
sion from the region of inelastic scattering and A is the 
absorbed part which is given by A = 1 — T co h — R co h- 
This model has been used by several other authors as 
well to simulate inelastic scattering. But, this model has 
a problem, in the presence of magnetic flux it is shown to 
violate the Onsager's two terminal symmetry relationsO, 
i.e., T tot ($) ^ Ttot(—$), where $ is the magnetic flux, 
and hence is not suitable in modeling inelastic scattering 
whereas Buttiker's model preserves the Onsager's sym- 
metry relations in the presence of magnetic flux. How- 
ever, the Buttiker model suffers from a shortcoming, in 
that it describes only localized dephasing (at the point 
contact between system and third lead) instead of de- 
phasing that occurs throughout the system in realistic 
situations. 

Brouwer and Beenakker! have removed the shortcom- 
ings of both the above stated models by mapping the 
three probe Buttiker's method into a two terminal ge- 
ometry, this is done by eliminating the transmission co- 
efficients which explicitly depend on the third probe by 
means of unitarity of the S-matrix. They have considered 
a three terminal geometry in which one of the probes is 
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FIG. 1: Aharonov 
three probe model. 



Bohm ring geometry. Inset shows the 



used as a voltage probe in absence of magnetic flux (see 
inset of Fig. 1). A current 1 = 1% = —I 2 flows from source 
to drain. In this model, a fictitious third lead connects 
the ring to a reservoir at chemical potential /13 in such a 
way that no current is drawn (I3 = 0). The 3X3 S- matrix 
of the entire system can be written as- 



S 



rn h2 ii3 

^21 r 22 t 2 3 
^31 ^32 ?~33 



Application of the relations! 



Hg],p = 1, 2,3 and G pq = {2e 2 /h)T pq 

yields the (dimensionless) two probe conductance G 
h 1 

2e 2 fli—H2 1 



G = T 21 



T23T31 
T31 + T- 



32 



(1) 



and on elimination of the transmission coefficients 
which involve the voltage probe using the unitarity of 
the S - Matrix leads toB 



G = T 21 



(1-R n -T 21 )(l-R 2 



T 21 ) 



1 — i?n — T 2 



1 — R 22 — Ti: 



(2) 



Now all the above coefficients are built from the 2X2 
S matrix- 



S' 



ni t 12 

hi r 22 



which represents the S Matrix of the absorbing system. 
The first term in Eq. 2 represents the conductance con- 
tribution from the phase coherent part. The second term 
accounts for electrons that are re-injected from the phase 



breaking reservoir, thereby ensuring particle conservation 
in the voltage probe model. For further calculation one 
can use the coefficients of the S' matrix of an absorbing 
system where absorption takes place uniformly. 

There are different ways of introducing absorption in 
the system -(f) The well known imaginary potential (IP) 
model and (2)Stochastic absorption (SA). In the first case 
Hamiltonian becomes non-hermitian due to the presence 
of complex potential term —iVi, which leads to absorp- 
tion (non conservation of particle mirnber). This model 
suffers from some spurious featuresEffl. In the scattering 
case, in the vicinity of the absorber, the particle experi- 
ences a mismatch in the potential (being complex) and 
therefore it tries to avoid this region by enhanced back 
reflection. Thus the IP model plays a dual role of an 
absorber as well as a reflector, i.e., absorption without 
reflection is not possible. Naively one expects the ab- 
sorption to increase monotonically as a function of Vi. 
However, the observed behavior is non-monotonic. At 
first absorption increases and after exhibiting a maximum 
decreases to zero as Vi — > 00. The absorber in this limit 
acts as a perfect reflector. During each scattering event 
an electron picks up an additional scattering phase shift 
due to Vi which along with multiple reflections leads to 
spurious scatterings (additional resonances )0tl in the sys- 
tem. As a result of the afore mentioned limitations there 
is need of a model wherein such spurious scatterings are 
absent and in the limit of large absorption does not cor- 
respond to a perfect reflector. Fortunately such a model 
exists, and this is the model of stochastic absorption!!!]. 
This model does not suffer from the drawbacks of the IP 
model as we will show towards the end of our work. 

The SA model is not new it has-aarlier been dealt with 
in the context of ID localization!!!. Stochastic absorp- 
tion in the ring is inserted by the factor e~ Q/l (or e _Q ' 2 ) 
in the comnlex free propagator amplitudes, every time 
we traversed the upper (or lower) arms of the ring (see 
Fig. I). We have calculated the relevant transmission 
and reflection coefficients by using the S matrix method 
along with the quantum wave guide theory for a single 
channel case. In this model, average absorption per unit 
length is given by 2a. With this method we show that 
the calculated conductance (G) in Eq. 2 is symmetric un- 
der the flux reversal as required. The visibility of the AB 
oscillations rapidly decay as a function of a, indicating 
dephasing. Hence forth we refer to a as an incoherence 
parameter. Increasing a corresponds to increasing de- 
phasing processes in the system or increase in tempera- 
ture. 

In Fig. 1, the length of the upper arm is l\ and that of 
lower arm is l 2 . The total circumference of the loop is L = 
l\ + l 2 . The loop is connected to two current leads. The 
couplers (triangles) in Fig. I which connects the leads and 
the loop are described by a scattering matrix S. The S 
matrix for the left coupler yields the amplitudes 0\ = 
(a',/3J,7^) emanating from the coupler in terms of the 
incident waves I\ = (a,/3i,7i), and for the right coupler 
yields the amplitudes 2 — (8' , /3' 2 , 72) emanating from 
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the coupler in terms of the incident waves I 2 = (S, fa, 72) 
The S-matrix for either of the couplerald is given fay- 



s' = 



1.0 




0.8 = 



with a = - 2e) - 1) and b = - 2e) + 1). 

e plays the role of a coupling parameter. The maximum 
coupling between reservoir and loop is e = i, and for 
e = 0, the coupler completely disconnects the loop from 
the reservoir. As mentioned above incoherence is taken 
into account by introducing an attenuation constant per 
unit length as mentioned before. 

The waves incident into the branches of the loop are 
related by the S Matrices for upper branch by- 
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FIG. 2: Conductance (G)for lengths h/l = 0.45, h/l = 0.55 
for coupling parameter e = 0.44 (waveguide coupling). Fermi 
wave-vector kl = 5.0. The inset shows the weak coupling case 
e = 0.10 for same physical parameters. 



These S matrices of course are not unitary S(a)S(ay 7^ 1 
but they obey the relation S(a)S(—a)^ — 1. The same 
relation is also obeyed by the S Matrix of the system in 
presence of imaginary potential. Here kl\ and kli are the 
phase increments of the wave function in absence of flux. 
^ and ^ are the phase shifts due to flux in the upper 
and lower branches. Clearly ^7- + ^ = where $ is 

the flux piercing the loop and $0 is the flux quantum^. 
The transmission and reflection coefficients in Eq. 2 are^ 
given as follows- T 21 = |^| 2 , R n = |^| 2 , i? 22 - |^| 2 ,;g 

T12 = |^-| 2 wherein 8' , S, a', a are as depicted in Fig. 1. ^ 

After calculating the required reflection and transmis- 
sion coefficients we graphically represent our results in 
the following figures. Throughout the discussion the 
physical parameters we give are in dimensionless units. 
We see that the coherent transmission T 2 i is not symmet- 
ric under the flux reversal however proper re-injection of 
carriers by Eq. 2 for the total conductance G plotted in 
Fig. 2 shows that the Onsager's symmetry relations are 
restored, i.e., G is symmetric under flux reversal. In the 
inset of Fig. 2 we have plotted G versus flux for the same 
physical parameters except we have considered the weak 
coupling case (e = 0.10). G is symmetric in flux and 
features are sharp as expected because of lifetime broad- 
ening of the energy levels is small in this case compared 
to the case for which e = 0.44. 

In Fig. 3 we plot visibility (V) as a function of inco- 
herence parameter a. Visibility is of course defined as- 
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The plot shows that with increase in the value of the 
parameter a the visibility exponentially falls off, reaching 



FIG. 3: Visibility for the same physical parameters as in fig- 
ure 2, coupling parameter e = 0.44. In the inset the harmonics 
have been plotted for the same physical and coupling param- 
eters. 



a point where it becomes zero corresponding to the disap- 
pearance of quantum interference effects. Iia-the inset of 
Fig. 3 we have plotted the first few Fourierliia harmonics 
Oj (wherein i = 1 to 4) of G(<&) as a function of a. The 
harmonics exponentially fall off with increasing a with 
exponent increasing as we go from 1 to 4. The nth order 
harmonic corresponds to the contribution from electronic 
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paths which encircle the flux n times. The harmonics 
can sometimes show non monotonic behavior depending 
on the physical parameters, however, the visibility is a 
monotonic function of the incoherence parameter a. 
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FIG. 4: Comparison of two models of dephasing. Herein the 
conductance (absorption in the inset) for the two cases SA 
and IP have been plotted. The parameters used are same as 
in Fig. 3. 

Now for the comparison between the two models, i.e., 
IP and SA, in Fig. 4 we plot the total conductance (G) 
for both the cases. In the case of imaginary potential 
we see that the total conductance (G) continuously de- 
creases and as Vi — ► oo, G — * (this limit is not shown 
in the graph as it goes beyond the scale used), while in- 



case of stochastic absorption as a — > oo, G goes to a 
constant value which depends on the Fermi-energy and 
other system parameters. Thus we see that the modeling 
of dephasing by SA is indeed justified as we want the AB 
oscillations to die out and not that the conductance itself 
should vanish, and this is where SA scores over the IP 
model. In the figure inset we have depicted the behavior 
of total absorption in the system A for both the cases 
for the same physical parameters. As a — > oo there is 
a finite absorption in the system as the electron prop- 
agates in the medium in this limit whereas absorption 
in the imaginary potential model is non-monotonic and 
in the limit Vi — > oo absorption vanishes. This is due 
to the fact that in the IP model as Vi — > oo absorber 
acts as a perfect reflector, there is no absorption in the 
medium as the particles do not enter the medium (and 
hence G = 0) obviously which is an unrealistic situation 
for real systems. 

In conclusion, we have shown that G(<&) is symmetric 
under flux reversal in the presence of incoherent scatter- 
ing represented by the incoherence parameter a. For this 
we have used the procedure of Brouwer and Beenakker in 
the presence of magnetic flux and simulated absorption 
(dephasing) using the method of stochastic absorption. 
We have used this method of SA to study the behav- 
iors of the various quantum phenomena, e.g., transport 
across resonant tunneling-systems, current magnification 
effect in mesoscopic ringsj!3 delay and sojourn timesE-3 in 
mesoscopic systems in the presence of incoherence. These 
results will be published elsewhere. 
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